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We study the condensate dynamics of the so-called entangled Bose-Einstein condensation (EBEC), 
which is the ground state of a mixture of two species of pseudospin-i atoms with interspecies spin- 
, exchange scattering in certain parameter regimes. EBEC leads to four inter-dependent superfluid 

components, each corresponding to the orbital wave function associated with a spin component of 
{^jfj' a species. The four superflows have various counter-relations, and altogether lead to a conserved 

r-^ , total supercurrent and a conserved total spin supercurrent. In the homogenous case, we also obtain 

the elementary excitations due to variations of the single-particle orbital wave functions, by exactly 
solving the generalized time-dependent Bogoliubov equations. There are three gapless Bogoliubov 
OO ', modes and one Klein-Gordon-like gapped mode. The origin of these excitations are also discussed 

from the perspective of spontaneous breaking of the symmetries possessed by the system. 
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I. INTRODUCTION 



PACS numbers: 03.75.Mn, 03.75.Gg 

' 

Bose-Einstein condensation (BEC) amplifies quantum mechanical behavior of individual particles into macroscopic 
quantum phases. In the simplest case, BEC occurs in a single-particle spatial or orbital state. Quantum features 
are more pronounced when there are additional degrees of freedom. For example, BEC occurring in a superposed 
G ' single particle state leads to Josephson effect [ij, while BEC of atoms with spin degree of freedom leads to a spinor 
condensate For spin = 1, in a mean field state 0,3; which is exact for a gas with ferromagnetic spin exchange, 
^ BEC occurs in a single-particle superposition of the thee hyperfine states. On the other hand, in the exact ground state 
O , for a spin-1 gas with antiferromagnetic spin exchange [4-7], BEC occurs in a superposition of two-particle state. As a 
further development in this perspective, the so-called entangled BEC (EBEC), i.e. BEC occurring in an interspecies 
I entangled two-particle state, amplifies entanglement of individual distinguishable particles into a macroscopic phase 
fS| • of a many-particle systems EBEC was found to be the ground state of a mixture of two species of pseudospin-i 
^ ' atoms in a considerable parameter regime (ol [l0|. It is entirely different from the two-component BEC, which occurs 
. in a gas of two kinds of atoms distinguished by only one degree of freedom 0, [ll| , whose ground state is simply a 
2^ ' direct product of the states of the two kinds of atoms, each separately undergoing BEC. Just like the simplest BEC 
may be a source of coherent atoms, EBEC could be a source of entangled atom pairs. 

EBEC, in the case where the total number of atoms of each species is equal to TV, refers to the many-body ground 
0\ •■ state 
O 

a^: ,^ . 1 

^/N + lN\' 



where and are, respectively, the annihilation operators of the two species a and b for pseudospin a (cr ='I",J,). 
In |Go), BEC occurs in a two-particle state of a maximally entangled interspecies pair 

V{^a,rb) = -^[0at(ra)| t)a0&i(rb)| i)b - <Pai{ra)\ i)a<Pbt{n)\ t)b], (2) 

and is thus called EBEC or BEC with an entangled order parameter ri{ra,Vb)- Here (paa (a = a, 6, cr —t,i) is the 
single-particle orbital wave function for each spin state of an atom of each species. In the most general case, spin 
dependence of the potential Uaa and the scattering lengths, which determine the orbital wave functions, lead to a 
kind of spin-orbit coupling. Consequently, ^ (f'ai: hence spin and orbital parts in 77(ra,rf,) cannot be factorized 
as an orbital part and a spin part, i.e. there exists spin-orbit "entanglement" in 77. If, however, (jja-f — (j)a^ — 4>on as 
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in the case where they are dominantly determined by the spin-independent part of the Hamihonian, then ?7(ra,r;,) = 
4>a{^a)4'b{i^b)^{\ t)a| i) b — I i)a| t)b), i-G- Spin and Orbit become disentangled, as in the usual consideration in SU(2) 
model spin-1 model, consequently, the entanglement between a-atom and 6-atom becomes entirely spin entanglement. 

A state analogous to \Go) appears in a SU(2) symmetric model of a single species of pseudospin-i atoms [1, 0], 
with the role of the two different species played by the single-particle orbital ground state and the first excited state 
of a single species of atoms. But there are also differences: (i) In the SU(2) model, the number of atoms in the two 
single-particle orbital states are not conserved, but are fixed by a measurement or controlled by a microcanonical 
distribution, while in EBEC the number of atoms in the two species are strictly conserved, (ii) In the SU(2) model, 
in which there is only a single species, the occupation of two orbital modes, rather than a single mode, is due to 
the constraint of spin conservation in the cooling process, which could be compromised as the collision rate during 
the evaporative cooling might depend on hyperfine spins, and there could be atom loss. In EBEC, in contrast, two 
distinguishable species can have small total spin, as distinguishable atoms are not subject to the constraint of Bose 
symmetry, consequently all the atoms of each species can occupy the lowest orbital modes, (iii) In the SU(2) model, 
the analog of |Go) is not the true ground state of the system. In EBEC, \Gq) is the true ground state, (iv) In the 
SU(2) model, the orbital identity of atoms, unlike different species, and the correlation between identical particles in 
the two orbital modes, unlike the entanglement between atoms of different species, are lost after the atoms are taken 
out of the trap. 

An important open issue about EBEC is its condensate dynamics, i.e. superfluid behavior determined by the orbital 
wave functions and elementary excitations due to the fluctuations of the orbital wave functions. In this article, under 
the presumption that EBEC exists, we study the condensate dynamics of EBEC, based on a generalized version of the 
time-dependent Gross-Pitaevskii (CP) equations. EBEC leads to four inter-dependent components of the superfluid 
with a few counter-relations between each other. We also study the elementary excitations by exactly solving a set of 
generalized version of the Bogoliubov equations, as well as from the perspective of symmetry breaking. 



II. HAMILTONIAN AND THE GENERALIZED GROSS-PITAEVSKII EQUATIONS 

Consider a dilute gas of two species of atoms in a trap. Each atom has an internal degree of freedom represented 
as a pseudospin-i. The field theoretic Hamiltonian is 

^ = E / d'r^Lhaai^ca + ^ ^ dll'^ j rf'^V'L V-L' ^a. + Uab. (3) 

Oca otaa' 

with 

^ab = E 5^"'' / d\^l,^l,^Pba'^aa + ge J Vfe^bt^ai + ^l^^ltMat), (4) 

where a ^ a,b represents the two species, a =tj-l' represents the two basis states of the pseudospin- ^ , haa = 
— ?;,^V^/2TOct + Uaa is the single particle Hamiltonian, ipaa is the field operator of species a. The coefficients g's are 
shorthands for g'^Sia-iai = {'2,Trh'^^i^ala3a4/ tJ^af}), where ^.jj^alasa^ is the scattering length for the scattering in which 
an Q;-atom flips from (74 to cr-i while an /5-atom flips from (T3 to fJ2, /^q/^ — Tfia 'mi3/{ma + mp) is the effective mass. 

For scattering lengths, we use the shorthands S^af^ = S^oTX, = 2^**1 = '^dZal ior a ^ a 1], ^J"^) = Ci°'a'a'a^ 

^{ab) ^ ^(ab)^ ^ ^(ab)^ ^ Correspondingly g'fj^ = 2Trh^ ^"J^ / ^lap , that is, giT^ = giZZl, glT'' = 2gi%%i = ^giZZl 

J. I - (ab) tab) (ab) (ab) fTTJ 

for CT ^ CT, g^^, = g„„,„,„, ge = fftw = H^- 

Under the single orbital mode approximation for each species, for each atom of species a {a — a, b) and pseudospin 
a (fj =t,i), only the single-particle spatial ground state 4>aa{j) is occupied. Therefore ipa = Q.a4>aa, where 
is the annihilation operator corresponding to the single-particle orbital wave function (j)aa- Then the many-body 
Hamiltonian can be simplified as 

H - ^ fa^aNa^a + ^ E ^l^fN^aNc^a' + ^ K^f} NaaNba' + A'elajai^^t + al^t^^i): (5) 
a, a a.aa' aa' 

where N^a — ci'a'^a- The total number of atoms of each species — Naf + Nai is conserved. The coefficients K^s 
are shorthands for 

KifaLa, ^ a^afXa.a, J C i^Wf^a, (r)0/3.3 (r)</>a.. (r)dV, (6) 
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that is, Kt^'^ ^ KiZZl, kIT^ ^ 2Kl%t = ^K^Zl for a ^ a, Kif) ^ K^^X'.^ ^ ^jl^i = <?t 0- 

faa = <^aa " i^ia^V^, where e^c = / (j)*a„haa4'aad^r IS the smgle-particle energy. 

When Na = Nb, the ground state is exactly \Go), under the following conditions, which ensure consistency of 
simplifying the Hamiltonian to that of isotropic Heisenberg coupling between two giant spins representing the two 
species: (i) Ua^{r) = Uai{r) — Ua{r). (ii) The intraspecies scattering lengths satisfy = ^aZ"^ — ^J^- (iu) 

The interspecies scattering lengths satisfy the relations ^j"''^ = ^j'^''^ , denoted as ^i"'''^ , and ^^"^^ = , denoted as 

^d"''^ , where the subscripts "s" and "d" stand for "same" and "different" , respectively, (iv) ^e"''^ = Ci"''^ — Cd"''^ • Under 
these conditions, fia-f — fJ-ai, (^at — 4'ai in the many-body ground state. The four wave functions (jiaa^s^ of which 
the condensate wave function 77 is built on, satisfy four generalized GP equations. It has also been shown that in a 
considerable parameter regime, the ground state approaches |Go). The conditions (i) and (ii) are also satisfied in the 
SU(2) symmetric model of a single species of pseudospin-i atoms. 

In the thermodynamic limit, the energetic advantage of EBEC is lost to the stability of simple BEG, like other 
fragmented BEG. Hence EBEG should be realized in a mesoscopic scale with a finite number of atoms and a finite 
volume. We can obtain upper bounds on the number of atoms as the following. According to previous discussions [lo| . 
the effect of symmetry-breaking perturbation, which causes the ground state to deviate from EBEG, tends to diminish 
when the volume remains finite while A < 2Ke , where A is the energy gap of the perturbed Hamiltonian, and is 
given by A = ^/W^, where d ee \ J,^-2K^^Ca~Cb\N^ + \Bb- Ba\N , N = {Na + Nb)/2, where = K^^^^+K^^^^- 

E^[K\f^+K^^^ - ifjf ) - ifg^)), with = - KZ, a = -HKj;^") - if}^") - i^l;-)), {a = a, h). Under 

the conditions given in the above paragraph, A 0, hence EBEG is indeed the ground state. With deviation from 
these conditions, we may use the requir ement A < 2Ke, i.e. d < K^, to derive a constraint on N for the occurrence 
of EBEG, which turns out to be iV < {y^\Bb - Ba\'^ + K^\J^ - 2K^ -Ca- Cb\ - \Bb - Ba\)/{2\J, - 2K^ -C^- Cb\). 
Furthermore, a necessary condition for EBEG to occur is that fc^T < Ei — Eq, where Ei and Eq are the energy of the 
first excited and the EBEG ground states. Ei — Eq 2Ke- According to dH), Ke — ge J ^*-|.(r)</)^j^(r)0b^(r)0a4,(i")rf^?'- 
To make a rough estimation using uniform wave functions, we have = ge/^, where fi is the volume of the gas. Hence 
we should have T < 2ge/kBfl. Alternatively, for atoms in a trap, we may roughly assume (j^aa ~ i^jw)^^^ exp{— "^^^ ), 
where it is assumed that the atoms of both species have equal mass, and that the trap is isotropic with frequency to. 
Then « gf,^mijj / {2ttK). Hence we have T < (5e/fcs)\/2ma;/(7r?i). BEG transition temperature can be roughly 
estimated using one component of the gas with particle number N/2 [l^. Hence Tc ~ 3.31?i^(A^/57)^/^/m for a 
uniform gas or Tc ~ 0.94?iajiV^/^ for a trapped gas. Moreover, combining the estimation of Tc with the result derived 
from ksT < Ei — i?o, we obtain a further constraint on TV, namely N < (2(7eTO)^^^/[(3.31?i^fcs)^/^f2^/^] for a uniform 
gas, or TV < {2gefm^/'^ /[k%{2Trujf/'^h^^'^] for a trapped gas. 

EBEG might be experimentally realizable by using two species of spin-1 alkali atoms in an optical trap with hyperfine 
states constrained in a two-dimensional subspace of | f) = |^ = 2,mF = 2) and \ ]^) = \F — l,mF = 1) [Ifl- In 
order that the spin-exchange scattering is energetically guaranteed, we should have — = £{,-(- — ef,j, [9|, where 
Eqct (a — a, b) is the single particle energy. However, for two different species of alkali atoms, the hyperfine splitting is 
different. A method of overcoming this difficulty is to apply a magnetic field such that ea-[ — (-ai, which is now the sum 
of the hyperfine splitting Aa plus the difference in Zeeman shift of the two hyperfine states {2ga^F=2 ^ ga.F=i)fJ'BB, is 
the same for the two species. For an alkali atom, ga p = [F{F + 1) + J{J + 1) — Ia{Ia + ^)]/F{F + l)j where J — 1/2, 
la is the nuclear spin of species a. Gonsider species a to be ®''Rb while species b to be ^^Rb. Then Aa = 6835AIHz, 
Ab = 3036MHz, la = 3/2, and h = 5/2 4]. It can be estimated that B = 0.325T. 

If EBEG exists, then what about its physical properties? Here we focus on its orbital dynamics, which is determined 
by the generalized time-dependent GP equation 

ot Iiria o 3 

where a ^ a represents the species other than species a, and a ^ represents the pseudospin opposite to g. 

These time-dependent GP equations can be obtained from the static GP equations ^], by replacing /Xq,^ as ihd/dt, 
and can also be justified by using the action principle 

5 [ Ldt = 0, (8) 
Jti 
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where 

f dr{{ih/2) Y,{'PladKaldt - ^^,d<l>lJdt) - (Go|H|Go)} (9) 

is the Lagrangian functional for jGo). The static GP-hke equations have been obtained minimization of the energy 
functional under |Go). This self-consistent determination of the equation of motion of the orbital wave function actually 
underlies the derivation of the simplest GP equation, for which the many-body ground state is {\ / y/N\){a^)^ \Q) . This 
methodology has also been used, e.g., by Ashhab and Leggett in studying the SU(2) model 0. 

In studying the orbital dynamics of the simplest BEC, it is assumed that the system remains in BEG though the 
condensate wave function is time-dependent. Likewise, in the present case of EBEC, the generalized time-dependent 
GP equations ([7]) presumes that the system remains in the many-body ground states in the form of |Go), though 
the corresponding orbital wave functions 4>a„ is time dependent. In other words, the system remains as BEG though 
the orbital wave functions are time dependent. It is this kind of dynamics that corresponds to superfluidity and is 
considered here. 



III. SUPERCURRENTS AND SPIN SUPERCURRENTS 



The number density of species a with pseudospin a is 

rw = {N/2)cpl„4>o,„, (10) 

while the supercurrent is 

h N , 



2nu 2 



where Vqct is the superfluid velocity of species a of pseudospin a. 
From the generalized time-dependent GP equation ([7]), we obtain 



9n (r,t) ^^ ^^^^^ 
ot 

with 

Saa- = -Sba -Sag- 5*60- ^ ^ ^^"^ Im{(j)lg(j)b„(j)as(l>*aa)^ (13) 

which is due to interspecies spin-exchange and acts as a source. Thus the supercurrent is not conserved individually in 
each pseudospin component of each species. Equation (|13p indicates the counter relations between the two pseudospin 
components of a same species, as well as those between two components with a same pseudospin and of two different 
species. 

The total supercurrent of each species is conserved: 

^^""^ + ""^^ V.(J., + J„,) = 0. (14) 

So is the total supercurrent for each pseudospin, 

d{naa + riba) 

dt 

Of course, the total supercurrent of the four components is conserved: 

d 



V • (Ja^+Jb,) =0. (15) 



+ V-^J„<. = 0. (16) 



dt 

a.{T a, a 

Furthermore, we can also define spin density and spin supercurrent for each species, 



= m_F(a,t)nat + WF(Q;,i)na4,, (17) 
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Ja='n^F{ct,t)^at+mF{a,i)3al, (18) 

where m_F(a,(T) denotes the hyperfine z component represented by pseudospin a for species a. 
The spin supercurrent of each species is not conserved: 

^ + V • = [{niFia, t) - mpia, l)]S^^. (19) 

If niFia, cr) = triFib, a), then the total spin supercurrent of the two species is conserved, 

^^^^^^ + V-(J: + J,^) = 0. (20) 
Thus we have a conserved total supercurrent as well as a conserved total spin supercurrent. 

IV. HYDRODYNAMICS 

In parallel to some discussions in Ref. here we derive some hydrodynamic equations from the generalized GP 
equations. Our discussion is restricted to zero temperature, hence does not require thermodynamic limit and also 
applies to a finite system in a trap. 

With (paa = /atre**"", whcrc the phase is not singular, Eq. ^ yields 

^ = -—V • (/LV$a.) - + ^^9- f^^f^^U^U^ sin($g^ + - - *aff), (21) 
at TOq in 

which is just the continuity equation expressed in terms of amplitudes and phases, and 

(A^-1) (gg) 2 I (ga)^2 2iV + 1 2 (22) 

- cos($a^ + - $as - $a^). 



/a 



According to ([TT|) . 



hence the gradient of p2l) becomes 



/a<T = — V$„^, (23) 
m 



^ + ^"^"""L)' (24) 

where /-/ --_^!_V2/?5 + l 2(iV-l) (aa) (jV-1) (aa) (jV-1) (gg) (2W+1) (aa) ^(e) 

wneie /Iqo- — -Ima^naa y ii-acr ^ U aa ^ 3 i/CTCT "•gcr T 3 ijcrS H-aa^ 3 i/crcr "gcr 1 3 i/crCT "■gcr ^gcr 7 

where the third to sixth terms on right-hand side are mean-field interaction energies due to scattering without spin 
exchange, while 



^(e) _ -f 2 Irinrrnnan^ 



^ gej „ COs($a. + $ag - ^ag - $a.) (25) 



is due to spin exc hange scattering. S^V^^ = [{N + 2) /3]geC^ \/naanaanaa /naa) cos{<Paa + *a<T - ^as - ^acr) - 

[{N + 2)/{'in)]geyJna9naanas/naa s\n{^aa + ^as - $a<T " ^'aa) [wa (va^r " ^&s) " maivaa ~ ^as)], implying the 
inter-dependence between superfluid velocities of different components. 

It is interesting that with the spin-exchange term, Eq. (I22p can still be written as a generalized Josephson relation 

^ = (26) 

dt h SUaa ' 

where 

E = Ea+Eb+ Eab, (27) 



6 



with 



3 9ei(t>tl(f>bt(t>*at^al + ^ai^aT'^bt'?^^)}' 

In case is independent of position or time, this leads to time- independent GP equations. 
Like the usual condensate, Eq. p4)) can be expressed in terms of pressure, 



(29) 



VFaa-V(^) + V(- —V'y^) VUaa, (30) 



dt nianaa 2 rria 2may/n^ nio 

where pressure Paa is related to chemical potential through 

dPaa = Uaadn'a^, (31) 

where fi'^^ = jlaa — Uaa- ^J■'aa ^^"^ thus Paa are contributed by interspecies spin-exchange scattering. 

V. ELEMENTARY EXCITATIONS 

The variation 64>aai away from in the ground state, satisfies the equation 

H 9cTa Was\ <30Q^ H ^ '(pQff'paffO^Qff H gijS 9aS(Paad<PaS 

, (2iV+l) |2c, , (2jV+l) (aa),* , , (2iV + 1) (.a), , 

iN + 2) {N + 2) {N + 2) 

ge(Paa<PasO(paa 9e(Paa'PaaO(paa ^ 9e(Paa'PaaO(paS , 

where (j)aa- refers to orbital wave functions in the ground state. 

For simplicity, here we only consider a homogeneous system, i.e. Uaa — 0- It is reasonable to suppose that 
a trapping potential Uaa would not change the essential physics. Furthermore, we assume the above-mentioned 
conditions giZ"^ = g^"^ = 9a, 9aa^ = gs, and g^^^^ = gd- Then the chemical potentials ^lai; = Hai, which is equal to 

N-l 2N+1 N + 2 

Thus single particle wave function in the ground state is 0Qcr = e^'^^'/'^/v^, which is then substituted to Eq. (|32|) . 

Even though c/ja-t = 4>a^i their variations should still be considered respectively. Since we are considering a uniform 
system, we may set 

K.(g)e'(''---"*) - t;L(g)e-('i--~-*)], (34) 



Then Eq. (|32]) yields 



w 2(iV-l) 

UJUaa = WUaa gaVaa 

,N-l N + 2 N-l 



I^^-l A^?2 (35) 



2iV-hl , 27V+1 N + 2 ^ 
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and 



where 



2(iV- 1) 

^Vaa = rrr 5qMqct - WVaa 

TV- 1 r 1 iV + 2 ^ 

+^-- + [-^- + ;^-]- (36) 

5(Ar_i) iV-1 27V+1 
^ = + 2^^ + ^17—^" + ^-9^ + ^^^^ 

h^q^ 2(N-l) TV + 2 

ZlTLa oil 6\l 

where the second equahty is a consequence of ([33]). 

Therefore we have eight coupled equations, which can be written as a matrix 

AV = wU, (39) 

where 

U = {Ua\, Wat, VaX, Ubf, W&f "H' "h)^' (40) 



the matrix elements of A can be read from Eqs.([35|) and 

We obtain four pairs of eigenvalues lu — ±.Eq'^ jh, (j = 1, 2, 3, 4), e\^^ being the energy of the elementary excitations. 
As usual, a positive eigenvalue cj corresponds to addition of a quasiparticle with momentum ?iq and removal of a 
quasiparticle with zero-momentum, while a negative eigenvalue a; corresponds to removal of a quasiparticle with 
momentum — ?iq and addition of a quasiparticle with zero- momentum. 

E^"* (j — 1,2,3,4) can be obtained exactly. The first and second excitations are Bogoliubov-like gapless modes, 
and can be written together as 

= ± [(^^)^ + (41) 

where j — 1,2 correspond to + and — in ± respectively, = h^q^/2ma, Xa = Eg [e^ + ^'"^^^^ ga], D = [{N + 2)ge — 
(N - l)gs - (27V + l)gd]'^. In terms of g, 

4^"^ - M(^+^)?^v+^(^+^)±[{(^-^)nv+^(^-^)r+,^^^]n*, (42) 

If ma = mh = m, then these two energy spectra reduce to 

4^)=i?°[i?° + r(^-)]}i, (43) 

where - h\y2m, r(^) = (iV - l){ga + gb)/n ± {{N - l)^{ga - gt)^ + AD/9y/yn. 

Prom p2)) . it can be seen that in the long wavelength limit 9 — >■ 0, the excitation Ejf'' (j — 1,2) becomes linear 

E(J) s^i^hq, (44) 

with the sound velocity 

Z\l ma mb 2s i TOa rrib yil'^mamb 

The existence of £'^^^ for q -> is subject to the condition {N — Vj^gagb > ^/9, or 50,9;, > [g^ ~ gs^ "^gdY ^ 1. 

In the short wavelength limit <? — >■ 00, 

ijW^,^, (46) 
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that is, they reduce to single particle excitations. 
The third and fourth excitations can be written as 

Ef=\Z,T{Zl-Y,)^, (47) 

where j = 3, 4 correspond to - and + in ± respectively, = e;^[e^/2 + (A^- l)ga/(3f^) + 2(A^ + 2).ge/(3fJ)] + e^[e^/2 + 
{N-l)gbl{:m) + 2{N + 2)gJ{in)]+R, where i? = 2(iV + 2)ge[(^-l)(ga+56-2ff.) + 2(2iV + l)5d + 2(7V + 2).9,]/(9r!2), 
Y, = e-e\{e-e\ + F) + 4e^(iV + 2)g,{[{N - l)gb + (TV + 2)g,]eym'') + 20/^ + 4e^(^ + 2)5e{[(iV - 1)5, + {N + 
2)r7e]e^/(9r!2) + 2C/3C!}, where F = 2[{N - l)gk + 2(iV + 2)ge]e^/(3n) + 2[{N - l)^, + 2{N + 2)ge\t\/{in) + 8(iV + 
2)2,g2/(902) + 4C, with C = {[(TV + 2)g, + (TV - l)ga][{N + 2)g, + {N - 1)5,] ^ [(TV - - (2TV + l)gd]2}/(9f]2). 

It is then clear that as g — >■ 0, — >• i?, Yq — >• 0, hence E^^ is gapless while e'^^ is gapped with the gap e'^^^ = \/2R, 
under the condition i? > 0, that is, {N -l)\ga+gb-2gs) + 2{2N +l)gd + 2{N +2)ge > 0, or ga+gb-2gs+4gd + '2ge > 
when TV > 1. 



In the long wavelength limit g — >■ 0, 



with the sound velocity 



5(3)^5, (48) 



^(3) _ ^2(TV + 2)5eC, 1 , 1 ,,1 



Eq"^^ being gapless is subject to the condition C/R > 0. When N ^ 1, s^^^ = [{ge + ga){9e+ Qb) — {Qs + 25d)^](l/?7ia + 
l/m6)/[3(sf„ + gb- 2g, + Ag^ + 2g^)n]. 
In the long wavelength limit g — )■ 0, 

Ef^ ^ yX^T^^^ (50) 
.... (51) 



2to, 



e// 



where 



A = V2R (52) 
is the energy gap at g = 0, similar to the energy due to the rest mass of a relativistic particle, 

2 _ N-l g^ gt^ 2(TV + 2) 5^ _ C 

~ 30 L„ + m6^+ 30 ^m^^ rub'^ r' ^ ^ 

is a constant similar to the square of the speed of light for a relativistic particle, 

A , , 

^eff = -2 (54) 



is the effective rest mass of the particle-like excitation near g = 0. Therefore, in the long wavelength limit, while the 
iree excitations Eq^\ E^^ and E^'' bet 
lein-Gordon particle, under the conditi 
In the short wavelength limit g ^ 00, 



three excitations E'"q\ E^^ and E^'' behave like massless particles, the fourth excitation E^^ behaves like a massive 
Klein-Gordon particle, under the condition i? > and (? > 0. 



Ef^^el, Ef^^e\ (55) 

if nia > nib- E^'' e^, E^^ if mb > nia- 

It is always the case that when g — )• 00, the excitations reduce to the free particle spectra and e^, each being 

double degenerate as there are two pseudospin states for each species. The effect of interspecies spin exchange is 
manifested as g — >■ 0, in both a nonzero sound velocity of the third excitation and the nonzero gap of the fourth 
excitation. 
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VI. NATURE OF SYMMETRY BREAKING 

From the point of view of symmetry breaking, Bogoliubov modes are Goldstone modes associated with U{1) gauge 
symmetry breaking. In our model, there are three conserved particle numbers Na-, iVf, and Na^+Ni,^ — Nai — Ni,^ = 23^, 
hence there are three C/(l) symmetries, the breaking of which gives rise to the three Bogoliubov modes. U(l) group 
is isomorphic to S0(2). Indeed, the U(l) symmetry generated by Na^ + Ntf — Nai — N^i is just the S0(2) symmetry 
generated by S^. 

It is not difficult to identify Eq^'' and Eq^"^ as the Goldstone modes associated Na and Nb, respectively. In fact, 

they reduce to the spectra of a-atoms and &-atoms if there is no interspecies scattering. Eq is associated with 

Na-^ + — Nai — Nfii- It is interesting that at the long- wavelength limit, Eq'^'' depends only on the spin-exchange 
scattering. 

The gapped mode -Eg^^ is due to the spin exchange between the two species. In fact, the gap vanishes if — 0. 
At the isotropic parameter point, the Hamiltonian can be rewritten as % = 2KeSa • Sb, where Sa is the total spin of 
a-species. Hence 

n = 2KeSaSb cos 9, (56) 

where Sa = Na/2, Sb = Nb/2, 9 is the angle between Sa and Sb- For the ground state, and 9 is uniquely tt. Around 
9 = 7:, the Hamiltonian is 

U, = 2K,SaSb[-l + {59)% (57) 

where 69 = 9 — tt. Therefore 69 is massive. This leads to the gapped mode. 

One can also consider this issue in terms of the phases ^aa- of the four components. In the absence of the interspecies 
spin-exchange part, "He, of the Hamiltonian, Naa- would be conserved. Then the spontaneous breaking of these four 
U{1) symmetries would give four phases $qo-'s, as well as four Goldstone modes, each corresponding to a combination 
of the four phases. With interspecies spin exchange, the Hamiltonian imposes an extra constraint on these four phases, 
as 

n = 2Ke{nai:nainb^nbiy/^ cos($at + *6i - - $bt)' (58) 
which fixes one combination of the four phases. By comparison of (I56p and (j58|) . it can be identified that 

9 = <i>a^ + $64 - <I>ai - (59) 

Around the minimum of T-le , 

= 2K,ina^nainb^nbiy/^[-l + (6^)% (60) 

where 5$ = ^af + ^bi — ^ai — ^bf — tt. Hence the mode corresponding to (5$ becomes gapped. This gapped mode 
corresponds to the source term Saa in Eq. (|12p . Therefore there are only three Goldstone modes remained. The 
gapped mode is due to the fluctuation oi 9 = $at + — ^ai — *&6f 

In the number conserved ground state \Gq), (a|a|) = = N/2, implying off-diagonal long-range order in 

each pseudospin component of each species. Equivalently, in the language of gauge symmetry breaking, U{1) gauge 
symmetry for each pseudospin component of each species is broken in the symmetry breaking ground state, i.e. 
(a,.,) = V^e^*°- ^ 0. 

On the other hand, EBEC means {a-fbi — a^b-^) ^ 0, as can be justified as follows. In the symmetry breaking ground 
state, {a^b^ - 04,6^) = (A^/2){exp[i($at + ^h)] - cxp[i($ai -I- $f,t)]}- ^s $at + - *a4, - $bt = ^ in the ground 
state, we have exp[j(<i>at + ^bi)] = — exp[z($ai + ^b^)]- Therefore {a-fbi — a^bf) = iVexp[i($a-|- + 'I'b^)]. Therefore, 
differing from the case of molecular BEC, we have (a^) 7^ 0, {b„) 7^ 0, as well as (a-[-&j^ — a^bf) 7^ 0, i.e. both atoms 
and the nonlocal interspecies entangled pairs are Bose-condensed. These field theoretical discussions are all consistent 
with the exact results in the particle-number conserved ground state above. 

VII. SUMMARY 

To summarize, we have studied condensate dynamics of EBEC, which exhibits peculiar superfluidity. EBEC 
leads to various counter relations among the four superfluid components. The total supercurrent and the total 
spin supercurrent are conserved. For the homogeneous case, we have also studied the elementary excitations due 
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to fluctuations of the four orbital wave functions around those in the ground state of the system. There are four 
modes, three of which are Bogohubov-hke while another is gapped. The three Bogoliubov modes correspond to the 
spontaneous breaking of the U{1) symmetries associated with three conserved particle numbers, while the gapped 
mode is associated with the spin exchange. Alternatively, the excitations can be understood as that the interspecies 
spin exchange gives mass to a certain combination of the four Bogoliubov modes corresponding to the four components 
of the system, hence there are only three gapless modes remaining. 

The emergence of massless elementary excitations of superfluids have been regarded as a paradigm showing how 
effective theory emerges from physics above the "Planckian" scale flS]. In this perspective, our model provides a way 
generating massive Klein- Gordon particles. 
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This work is supported by National Science Foundation of China (Grant No. 10674030), Shuguang Project (Grant 
No. 07S402) and Ministry of Science and Technology of China (Grant No. 2009CB929204). 
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